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This  paper  focuses  on  trajectory  control  of  the  6-DOF  body  fixed  reference  frame  lo¬ 
cated  on  a  very  flexible  aircraft.  The  6-DOF  equations  of  motion  of  a  reference  point  on  the 
aircraft  are  coupled  with  the  aeroelastic  equations  that  govern  the  geometrically  nonlinear 
structural  response  of  the  vehicle.  A  low-order  strain-based  nonlinear  structural  analy¬ 
sis  coupled  with  unsteady  finite-state  potential  flow  aerodynamics  form  the  basis  for  the 
aeroelastic  model.  The  nonlinear  beam  finite  element  structural  model  assumes  constant 
strain  over  an  element  in  extension,  twist,  and  in/out  of  plane  bending.  The  geometrically 
nonlinear  structural  formulation,  the  finite  state  aerodynamic  model,  and  the  nonlinear 
rigid  body  equations  together  provide  a  low-order  complete  nonlinear  aircraft  analysis 
tool.  Due  to  the  inherent  flexibility  of  the  aircraft  modeling,  the  low  order  structural  fre¬ 
quencies  are  of  the  same  order  as  the  rigid  body  modes.  This  creates  a  coupling  which 
cannot  be  separated  by  previous  control  schemes.  The  flexibility  must  be  accounted  for 
directly  in  the  controller  development.  To  accomplish  this  a  heuristic  approach  based  upon 
pilot  behavior  is  developed.  This  approach  separates  the  problem  into  two  parts:  a  fast 
inner-loop  and  a  slower  outer-loop.  Dominant  kinematic  nonlinearities  are  handled  in  the 
outer-loop  while  the  inner  loop  is  further  separated  into  a  lateral  and  longitudinal  motion. 
Control  of  the  inner-loop  lateral  motion  is  accomplished  using  a  standard  Linear  Quadratic 
Regulator.  For  the  longitudinal  motion  Dynamic  Inversion  is  utilized.  Differences  between 
the  desired  and  actual  trajectories  are  handled  in  the  nonlinear  outer-loop  using  tradi¬ 
tional  proportional-integral-derivative  design  guidelines.  The  closed  loop  time  integration 
is  accomplished  using  an  implicit  modified  Newmark  method.  Numerical  simulations  are 
presented  highlighting  the  strengths  and  weaknesses  of  the  method. 

I.  Introduction 

On  17  December  1903,  the  Wright  brothers  set  about  the  task  of  launching  an  airplane  into  the  sky. 
Their  multiple  successes  that  day  have  been  hailed  as  the  start  of  heavier  than  air  powered  flight.  One  of  the 
key  features  of  their  aircraft  was  the  use  of  wing  flexibility  for  roll  control.  Due  to  the  low  dynamic  pressure 
seen  on  that  flight  and  the  relatively  high  stiffness  to  mass  ratio  of  the  aircraft,  the  Wright’s  were  able  to 
develop  the  required  control  power  without  any  detrimental  aeroelastic  effects.  However  almost  a  100  years 
later  on  26  June  2003,  NASA’s  Helios  aircraft,1 

HP03-2  took  off  at  10:06am  local  time  from  the  Navys  Pacific  Missile  Range  Facility  (PMRF) 
located  on  the  island  of  Kauai,  Hawaii.  ...At  10:22am  and  10:24am,  the  aircraft  encountered 
turbulence  and  the  wing  dihedral  became  much  larger  than  normal  and  mild  pitch  oscillations 
began,  but  quickly  damped  out.  At  about  30  minutes  into  the  flight,  the  aircraft  encountered 
turbulence  and  morphed  into  an  unexpected,  persistent,  high  dihedral  configuration.  As  a  result 
of  the  persistent  high  dihedral,  the  aircraft  became  unstable  in  a  very  divergent  pitch  mode 
in  which  the  airspeed  excursions  from  the  nominal  flight  speed  about  doubled  every  cycle  of 
the  oscillation.  The  aircrafts  design  airspeed  was  subsequently  exceeded  and  the  resulting  high 
dynamic  pressures  caused  the  wing  leading  edge  secondary  structure  on  the  outer  wing  panels  to 
fail  and  the  solar  cells  and  skin  on  the  upper  surface  of  the  wing  to  rip  off.  The  aircraft  impacted 
the  ocean  within  the  confines  of  the  PMRF  test  range  and  was  destroyed.  . . .  The  root  causes  of 
the  mishap  include:  [a]  Lack  of  adequate  analysis  methods  led  to  an  inaccurate  risk  assessment 
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of  the  effects  of  configuration  changes  leading  to  an  inappropriate  decision  to  fly  an  aircraft 
configuration  highly  sensitive  to  disturbances  [,  and]  configuration  changes  to  the  aircraft,  driven 
by  programmatic  and  technological  constraints,  altered  the  aircraft  from  a  spanloader  to  a  highly 
point-loaded  mass  distribution  on  the  same  structure  significantly  reducing  design  robustness  and 
margins  of  safety. 

The  Helios  accident  highlighted  our  limited  understanding  and  limited  analytical  tools  necessary  for  de¬ 
signing  very  flexible  aircraft  which  have  and  potentially  exploit  aircraft  flexibility.  The  number  one  root 
cause/recomendation  from  NASA1  was  a 

[That]  more  advanced,  multidisciplinary  (structures,  aeroelastic,  aerodynamics,  atmospheric,  ma¬ 
terials,  propulsion,  controls,  etc)  time-domain  analysis  methods  appropriate  to  highly  flexible, 
morphing  vehicles  [be  developed]. 

Despite  the  lack  of  fundamental  understanding  on  the  behavior  of  these  vehicles,  recent  advances  in  air¬ 
borne  sensors  and  communication  packages  have  brought  the  need  for  high-altitude  long-endurance  (HALE) 
aircraft.  These  platforms  can  be  categorized  under  three  broad  missions,  supporting  either  the  military  or 
civilian  communities.  The  missions  include  airborne  Intelligence,  Surveillance,  and  Reconnaissance  (ISR), 
for  the  military,2  network  communication  nodes  for  the  military  and  civilian  usage,3  and  general  atmospheric 
research.3  Due  to  the  mission  requirements,  the  desired  vehicles  are  characterized  by  high-aspect-ratio  wings 
and  slender  fuselages,  resulting  in  very  flexible  vehicles.  Examples  of  mission  optimization  studies  for  this 
class  of  vehicle  can  be  found  in  Ref.  2,  where  it  is  shown  the  aircraft  are  required  to  have  a  fuel  fraction 
greater  than  66%.  This  results  in  a  very  small  structural  weight  fraction.  The  combination  of  high  aerody¬ 
namic  efficiency  and  low  structural  weight  fraction  yields  inherently  flexible  wings  and  nonlinear  structural 
and  flight  dynamics.  The  HALE  vehicle  will  then  be  susceptible  to  large  dynamic  wing  deformations  at  low 
frequencies,  presenting  a  direct  impact  into  the  flight  dynamic  characteristics  of  the  vehicle  and  controller 
design,  as  was  seen  in  the  Helios  flight  tests.1 

The  mission  of  the  HALE  aircraft  is  planned  to  be  unmanned  due  to  its  “dull,  dirty,  or  dangerous”4 
nature,  i.e., 

the  attributes  that  make  the  use  of  unmanned  preferable  to  manned  aircraft  . . .  [are]  in  the  case 
of  the  dull,  the  better  sustained  alertness  of  machines  over  that  of  humans  and,  for  the  dirty  and 
the  dangerous,  the  lower  political  and  human  cost  if  the  mission  is  lost,  and  greater  probability 
that  the  mission  will  be  successful.  Lower  downside  risk  and  higher  confidence  in  mission  success 
are  two  strong  motivators  for  continued  expansion  of  unmanned  aircraft  systems. 

For  all  the  reasons  stated,  a  better  understanding  of  flight  dynamics  and  control  of  these  vehicles  is 
required.  This  paper  presents  a  control  architecture  specifically  designed  for  the  very  flexible  aircraft. 

A.  Previous  Work 

In  1914  Lawerence  Sperry,  son  of  Elmer  Sperry,  demonstrated  his  father’s  autopilot  over  Paris,  by  standing 
up  in  the  cockpit  of  his  airplane  and  having  his  mechanic  walk  out  on  the  wing  to  create  an  external  dis¬ 
turbance.5  Sperry’s  invention  was  capable  of  maintaining  pitch,  roll,  and  heading  angles.6  This  remarkable 
demonstration  ushered  in  the  use  of  the  autopilot  for  a  variety  of  aircraft  and  aircraft  missions.  Sperry’s 
autopilot  relied  upon  linear  techniques  and  a  relatively  stiff  biplane  box  type  of  construction.  Since  then 
control  theory  has  evolved  significantly.  For  an  overview  of  aeroelastic  control,  Mukhopadhyay '  has  provided 
an  excellent  review  paper  of  both  analysis  and  control  of  aeroelastic  structures  over  the  past  100  years.  He 
provides  a  summary  of  elastic  theory,  unsteady  aerodynamics,  control  and  highlights  the  future  of  nonlin¬ 
ear  aeroelastic  control.  Recently  researchers  have  looked  at  applying  linear,  robust  control,  and  nonlinear 
techniques  for  aeroelastic  control. 

Modern  linear  control  techniques  such  as  Linear  Quadratic  Regulator  (LQR)  and  Linear  Quadratic  Gaus¬ 
sian  (LQG)  observers  have  been  applied  by  several  researchers.  Newman  and  Buttrill8  have  researched  the 
longitudinal  flight  dynamics,  longitudinal  linear  controller  design,  and  sensor  location  for  a  high  speed  trans¬ 
port.  Their  aircraft  modeling  shows  the  importance  of  including  elastic  modes  and  sensor  location  when 
designing  controllers  for  flexible  aircraft  and  the  difficulty  in  suppressing  aeroelastic  dynamic  effects  from 
the  rigid  body  response  of  a  flexible  supersonic  transport.  Tuzcu9  and  Meirovitch  and  Tuzcu10  have  coupled 
nonlinear  rigid  body  dynamics  with  linear  structural  dynamics.  The  formulation  treats  structural  dynamics 
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II.  Theoretical  Formulation 


There  are  two  fundamental  approaches  when  developing  a  controller.  The  first  is  to  treat  the  system  of 
interest  as  a  black  box  where  the  physics  of  the  box  are  not  known,  but  rather  linearized  transfer  functions 
can  be  fit  to  input/output  relationships.  The  second  approach  is  to  start  with  the  known  physics  of  the 
model  and  develop  the  controller.  The  later  is  used  here  and  hence  requires  a  basic  understanding  of  the 
underlying  governing  differential  equations. 

The  objective  of  the  controller  is  to  provide  closed  loop  reference  tracking  of  a  body  fixed  reference  frame, 
B,  at  point  O,  which  in  general  is  not  the  aircraft’s  center  of  mass,  Figure  1 ,  while  including  the  effects  of 
nonlinear  aeroelasticity.  The  B  reference  can  either  be  fixed  inside  a  rigid  section  of  the  fuselage  or  more 
generally  be  attached  to  a  node  of  an  elastic  fuselage.  When  the  B  frame  is  attached  to  a  node  of  an  elastic 
member,  the  y  axis  is  chosen  to  be  tangent  to  the  undeformed  longitudinal  axis  of  the  fuselage.  The  x  axis 
is  chosen  to  be  positive  out  the  right  wing,  such  that  for  an  undeformed  wings  level  aircraft  orientation, 
the  x  —  y  plane  is  parallel  to  inertial  frame  G.  The  2  axis  is  simply  the  cross  product  of  the  x  and  y  axes. 
The  tracking  will  consist  of  maintaining  desired  linear  and  angular  velocities  of  the  B  reference  frame.  The 


Figure  1.  Basic  body  reference  frame  and  vehicle  coordinates 

means  for  propagating  the  reference  frame,  B,  forward  in  time  is  done  by  deriving  and  integrating  a  series 
of  first-order  differential  equations  of  the  form 


x  =  f(x,u)  (1) 

where  x  represents  the  states  of  the  reference  frame,  B ,  and  u  represents  control  surface  and  external  inputs. 
Depending  on  the  fidelity  of  the  analysis,  these  first-order  differential  equations  vary  in  their  complexity 
from  simple  linear  time  invariant  to  nonlinear  time  varying  differential  equations.  For  the  classic  “rigid” 
body  analysis,27  the  first-order  differential  equations  take  the  form 

vb  —  fvD(vB,^B,C>PB,go,‘m,  Fext) 

UB  =  fu>D(VB,lB,C,PB,Mext)  (2) 

C  =  fduB,  0 
PB  =  /pd(C^b) 
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as  disturbance  inputs  to  the  rigid  body  dynamics.  Results  are  presented  for  light  weight  transport  aircraft 
using  standard  LQR  and  LQG  control  theory.  Pedro  and  Bigg11  have  also  applied  traditional  LQR  and  LQG 
techniques  to  a  flexible  aircraft.  They  developed  a  simulation  to  evaluate  pilot  ratings  for  an  aeroelastic  lon¬ 
gitudinal  aircraft  model.  Their  approach  looks  only  at  the  longitudinal  dynamics  with  a  Dry  den  gust  model. 
Cesnik  and  Ortega-Morales12  have  utilized  the  geometrically  exact  beam  formulation  coupled  with  finite 
state  aerodynamics  and  imbedded  actuation  to  study  and  control  flutter  of  a  HALE  type  cantilevered  wing. 
Their  work  focuses  on  the  use  of  LQG  observers  to  optimize  sensor  type,  sensor  placement,  and  actuation 
distribution.  They  have  shown  that  a  single  strain  gauge  sensor,  optimally  placed,  can  provide  the  sensing 
necessary  to  control  flutter  of  a  high  aspect  ratio  cantilevered  wing. 

Robust  techniques  using  H\,  H2,  #00 >  and  /z-synthesis  have  been  studied  by  several  researchers.  Chavez 
and  Schmidt13  utilized  linear  structural  modes  to  develop  longitudinal  controllers  based  upon  /z-synthesis. 
Their  work  focuses  on  developing  a  systematic  approach  where  robust  controllers  could  be  controlled  without 
complete  knowledge  of  the  flexible  system.  They  also  present  techniques  for  identifying  structural  modes  of 
interest  and  methods  of  incorporating  unsteady  aerodynamics  directly  into  the  robust  formulation  without 
the  necessity  of  building  an  aerodynamic  state  space  model.  Kron  et.  al.14  have  applied  a  reduced  order 
Hoo  controller  to  a  two-DOF  lateral  model  of  highly  flexible  transport  aircraft.  Their  work  focuses  on 
characterizing  model  uncertainty  when  reducing  the  structural  model.  The  uncertainty  is  then  handled 
systematically  through  the  development  of  a  two-DOF  H0 0  lateral  controller.  Finally  numerical  results  are 
presented  validating  their  method.  Li  and  Agarwal15  have  also  studied  the  use  of  a  reduced  order  structural 
modeling  to  develop  a  linear  robust  controller  using  H 2  and  techniques  for  a  high  speed  civil  transport 
aircraft.  They  have  shown,  for  that  class  of  aircraft,  the  model  reduction  provides  sufficient  information 
for  the  linearized  controller.  Goman  et.  al.16  have  conducted  a  parametric  study  of  various  H ^-based 
and  traditional  LQR/LQG  controllers  on  a  longitudinal  elastic  aircraft  model.  Their  study  shows  similar 
controller  designs  between  the  two  methods  and  improved  robustness  using  the  //oo-based  design. 

Dardenne  and  Ferreres17  have  designed  a  lateral  controller  for  highly  flexible  transport  aircraft.  The 
design  assumes  a  linear  time  invariant  (LTI)  plant  and  uses  a  linear  quadratic/programming  procedure  to 
design  the  controller  utilizing  frequency  domain  constraints.  Their  initial  results  have  looked  at  reducing  the 
wing  bending  effects  on  lateral  motion.  Patil18  and  Patil  and  Hodges19  have  also  used  minimization  routines 
in  designing  a  Static  Output  Feedback  (SOF)  controller  used  for  flutter  suppression  of  HALE  type  aircraft. 
The  resulting  controller  is  of  a  much  lower  order  than  LQR/LQG  type  controllers,  but  is  valid  only  for  a 
single  operating  condition.  The  authors  comment  that  due  to  the  simple  design  of  the  SOF  controller,  gain 
scheduling  of  different  SOF  controllers  might  be  easily  implemented. 

Calise  et.  al.20  have  studied  the  use  of  output  feedback  (dynamic  inversion)  mixed  with  a  neural  net¬ 
work  plant  perturbation  estimate  to  control  a  longitudinal  flexible  aircraft  model.  The  process  assumes  a 
non-minimum  phase  (NMP)  system  with  stable  zero  dynamics.  They  present  initial  results  showing  that 
their  method  could  suppress  structural  mode  interaction  for  a  simplified  longitudinal  aircraft  model.  Kr- 
ishnaswamy  and  Bugajski21  have  used  dynamic  inversion  for  studying  control  of  rockets  with  fuel  slosh.  A 
key  aspect  of  their  study  is  the  use  of  an  underactuated  system.  They  have  developed  an  observer  model 
for  estimating  the  fuel  slosh  dynamics  and  show  that  the  resulting  controller  can  control  the  pitch  dynamics 
of  booster  vehicle  while  damping  fuel  slosh  dynamics.  Finally,  Gregory22*25  has  applied  a  “novel”  Dynamic 
Inversion  control  technique  for  suppressing  longitudinal  motion  due  to  linear  aeroelastic  effects  of  high  speed 
transports.  The  variation  in  the  Dynamic  Inversion  is  the  inclusion  of  pre-filters  to  move  the  structural 
modes  further  into  the  left  half  plane.  Results  show  the  ability  to  significantly  improve  ride  control  and 
longitudinal  aircraft  handling  in  the  presence  of  longitudinal  structural  modes. 

While  most  of  the  aeroelastic  control  work  summarized  here  deals  with  different  classes  of  aeroelastic 
effects  than  the  very  flexible  HALE  type  aircraft,  the  control  techniques,  both  linear  and  nonlinear,  have 
potential  relevance  to  trajectory  control  of  very  flexible  aircraft. 

B.  Present  Work 

The  objective  of  this  paper  is  to  present  a  control  scheme  for  trajectory  tracking  of  very  flexible  aircraft.  It 
will  be  accomplished  through  a  coupled  6-DOF  vehicle  dynamics  with  a  modified  version  of  the  nonlinear 
strain-based  structural  formulation26  for  high-aspect-ratio  lifting  surfaces.  The  proposed  control  scheme 
is  based  loosely  on  human  pilot  operations.  Simulations  are  presented  both  symmetric  and  asymmetric 
maneuvering  and  different  aircraft  loadings. 
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where  the  B  reference  frame  linear  and  angular  velocity  variables  are  represented  by  vb  and  u Fext  and 
Mext  are  in  general  state  dependent  external  forces  and  moments;  m  is  the  aircraft  mass,  and  Ib  is  the 
aircraft’s  inertia  matrix  about  the  origin  of  the  B  reference  frame.  The  orientation  of  the  B  reference  frame 
is  accomplished  in  a  variety  of  ways  from  a  minimum  representation  using  three  non-orthogonal  Euler  angles, 
to  non-minimum  four  parameter  quaternion  representation,  to  a  nine  parameter  set  corresponding  to  the  nine 
components  of  the  set  of  unit  vectors  defining  the  triad  at  B.  Reference  28  provides  a  summary  of  different 
methods  used  in  the  aerospace  industry.  In  this  paper,  all  three  techniques  are  used  to  simplify  the  equations 
where  necessary.  In  Eq.  2,  £  is  the  vector  of  four  quaternion  elements  used  to  determine  the  orientation  of 
the  B  reference  frame  and  ps  is  the  inertial  position  of  the  B  reference  frame.  The  gravitational  field  effects 
are  represented  by  g0. 

A.  Summary  of  Governing  Differential  Equations 

The  rigid  body  formulation  has  three  key  assumptions  which  render  invalid  when  dealing  with  very  flexible 
vehicles:  1)  inertia  properties  are  constant  or  at  best  slowly  time  varying,  2)  the  coupling  inertial  force  due  to 
a  rotating  coordinate  frame  and  relative  velocity  of  flexible  members  is  negligible,  and  3)  external  forces  and 
moments,  Fext  and  Mext,  which  come  from  aerodynamic  loading,  are  based  upon  a  fixed  aircraft  geometry. 
In  the  rigid  body  case,  Eq.  2  presents  only  inertial  and  external  forces  and  moments.  For  the  flexible  aircraft 
a  set  of  elastic  EOM  is  also  introduced,  which  in  the  context  of  this  study  results  in 

Mq  +  Cq  +  Kq  =  R(q,  q,  A)  (3) 


«  1 

1 

«  1 

1  1 

*  I 

VB 

<?  = 

VB  | 

Q=\ 

VB  / 

.  eB  J 

1 

,  ub  j 

1 

,  J 

where  M  represents  generalized  mass  properties,  q  is  a  set  of  generalized  coordinates  containing  both  strain, 
6,  associated  with  the  flexible  vehicle,  and  the  inertial  position,  ps,  and  an  arbitrary  orientation  vector,  0b, 
of  the  B  reference  frame.  The  matrix  C  contains  both  structural  damping  and  nonlinear  terms  associated 
with  relative  position  and  velocity  terms  associated  with  a  rotating  coordinate  frame  (ujbx  vb,  etc),  K  is 
the  stiffness  matrix,  and  R(q,q,\)  represents  generalized  forces  (including  aerodynamic  forces)  which  are  a 
function  of  the  finite  state  inflow29,30  velocities,  A.  Coupling  of  the  rigid  body  and  flexible  dynamics  occurs 
through  the  dependency  of  M,  C,  and  R.  Typically  the  B  reference  frame  linear  and  angular  velocities  are 
represented  by 

<5> 

The  present  work  uses  a  constant  strain-based  formulation26,31  which  allows  for  airframe  nonlinear  geometric 
deformation  and  accounts  for  geometry-dependent  inertia  properties  of  the  aircraft. 

To  develop  the  nonlinear  governing  differential  equations  for  slender  elastic  structures,  a  systematic 
approach  is  used  where  the  rigid  body  and  elastic  EOM  are  developed  about  the  B  reference  frame.  The 
differential  equations  for  the  orientation  and  displacement  of  the  B  reference  frame  are  appended  based 
upon  a  four  state  quaternion  representation.  Unsteady  aerodynamic  modeling  is  included  and,  if  required, 
algebraic  equations  for  absolute  or  relative  constraints  are  appended  (example  of  which  is  a  joined  wing 
aircraft  as  shown  in  Figure  2,  where  relative  constraints  are  needed  at  the  joint  of  the  two  wings). 

The  derivation  of  the  EOM  is  based  upon  the  principle  of  virtual  work.  The  method  accounts  for  the 
virtual  work  associated  with  the  B  reference  system,  flexible  aircraft  slender  (beam)  structural  members, 
and  rigid  bodies  attached  to  the  flexible  structures.  The  virtual  work  of  a  beam  and  rigid  bodies  attached 
to  a  beam  are  initially  written  in  terms  of  dependent  displacement  vectors.  Then  the  kinematic  relationship 
between  beam  dependent  position  vectors  and  the  associated  strains  is  developed.  The  components  of  virtual 
work  are  summed  and  the  resulting  set  of  equations  are  transformed  from  a  set  of  dependent  position  vectors 
and  a  nonminimum  set  of  B  reference  frame  components  to  an  independent  set  of  strain  variables  and  body 
linear  and  angular  velocities. 

The  final  virtual  work  expression,  which  includes  both  B  reference  frame  and  flexible  body  contributions, 
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Figure  2.  Joined  wing  aircraft  concept 


is  written 


6W 


[  SeT  8bT  ]T 


Mff  Mfb 
Mbf  Mbb 


Kff  Kfb 
Kbf  Kbb 


where  the  generalized  mass  and  damping  matrices  are  given  by 


Cff  Cfb 
Cbf  Cbb 


Mff  =  J he  Mg  J he  Mfb  =  Jh€MGJhb 

Mbf  =  JhbMcJht  MBb  =  JhbMGJhb  +  MB 

Cff  =  J he  Mg  J he  +  Cg  Cfb  =  JheMcHhb  +  2J heMcHhi q 

Cbf  =  JhbMGJhe  Cbb  =  J^MG^hb  +  2JZbMGHh.p  4-  Cb 

Kff  —  Kg  Kfb  =  0 

K  bf  —  0  K  bb  —  0 


(6) 


(7) 


and  Mg ,  CG,  and  KG  are  the  assembled  flexible-element  generalized  mass,  damping,  and  stiffness  matrices, 
respectively.  The  matrices  Mb  and  Cb  represent  the  mass  and  damping  matrices  associated  with  the 
B  frame  rigid  element  portion  as  described  in  Refs.  32,33.  The  jacobian  matrices,  Jhe  and  Jhb,  provide 
relationships  between  flexible  position  and  orientation  vectors  and  the  independent  coordinates  of  strain, 
e,  and  B  reference  frame  linear  and  angular  velocities,  (3.  Additional  matrices,  Hhb  and  capture  the 

dynamics  of  a  rotating  coordinate  frame.  Complete  details  of  the  derivation  are  provided  in  Refs  32,33. 
Note  the  traditional  aircraft  rigid  body  EOM  can  be  recovered  from  the  flexible  EOM,  Eq.  (>,  by  assuming 
the  elastic  DOF  are  constant.  This  is  described  in  detail  in  Ref.  33.  The  resultant  force  vector  R  is 


R  = 


Rf 

Rb 


} 


Kff 

Kbf 


£ initial 


BgF 

B  g  B 


9B  + 


Bp, 

Bpo  j 


Fpt  + 


dM 

B™ 

Pb  J 


Mpt 


BLF 

dF 

adstD 


pdst  _j_ 


B&f 

B&. 


\ldat  + 


where  e initial  >s  an  initial  strain  vector,  gB  is  the  body-frame-  H  resolved  gravity  vector,  and  Fdst ,  Mdst ,  Fpt, 
and  A/pt,  are  body  resolved  distributed  and  point  forces  and  moments.  The  aerodynamic  forces  and  moments, 
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p°ero  ancj  faero ,  wjjjch  are  functions  of  control  surface  inputs,  u,  are  included  in  Fdst  and  Mdst.  The 
influence  matrices  are  derived  as31 

bLf  =  J£bf 
BLf  =  J$(Bm 

Bf  =  JT 

^VF  Jpt 

pM  _  jT 

aPF  -  J0t 

where  the  matrices  BF  and  Bm  are  constant  matri 
For  more  details  please  see  Ref  31.  For  simplicity  the  initial  strain,  e^mai  and  point  moments,  Alpt  will  be 
assumed  to  be  zero.  The  finite  strain  formulation29,30  for  aerodynamic  forces  and  moments  is  linear  in  the 
discrete  trailing  edge  surface  deflections.  These  assumptions  and  formulations  combined  with  a  point  force 
representing  a  simplistic  engine  model  allow  the  generalized  force  to  be  expressed  as 


BLo  =  Jp0BF 

(8) 

Bdlto  ~  JosBm 

(9) 

r  ^ 

II 

CQ 

(10) 

ja 

11 

(11) 

defined  by  an  elastic  element’s 

undeformed  length. 

p  pctero 

R  =  \Bg\  gB  +  [fl£t]  Faero  +  [BZ]  Maero  +  [BZ\ - -  r  dF1 


OUflap 

For  additional  details  see  Ref.  26.  From  the  Principle  of  Virtual  Work,  Eq.  6  yields 


Uflap  \Bp  ]  Uthrust 


Mff 

Mfb 
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4- 

Cff 

Cfb 

e 

4- 

Kff 

Afb 

6 

.{ R" ) 

Mbf 

Mbb 

.  P  . 

Cbf 

Cbb 

.  0 . 

Kbf 

Kbb 

b 

1  Rb  1 

(12) 


(13) 


This  set  of  equations  are  the  ones  given  in  a  compact  form  by  Eq.  3.  Note  that  M  =  M(c)  and  C  =  C( c,  c,/?), 
and  K  the  generalized  stiffness.  All  the  other  nonlinearities  are  contained  in  the  generalized  force,  R. 
When  the  EOM,  Eq.  13,  are  augmented  with  the  B  reference  frame  orientation,  position,  and  unsteady 
aerodynamics,  the  complete  set  of  governing  differential  equations  is 


(14) 


where  £  is  a  vector  of  four  quaternion  parameters  used  for  the  orientation  of  the  B  reference  frame,  £b  is 
the  the  inertial  position  vector  of  the  B  reference  frame  (Figure  1),  CBG  is  a  transformation  matrix  between 
a  B  reference  frame  vector  and  an  inertial  ((7)  vector,  A  is  a  set  of  unsteady  aerodynamic  inflow  velocities 
with  associated  differential  equation  matrices  F\  through  F3.  For  complete  details  please  see  Refs.  32,33. 


MfF €  — 

—MfbP  —  Off*  —  Cfb  ft  —  Kff t  4-  Rf 

MbbP  = 

—  Mbf £  —  Cbb  ft  —  Cbf *  +  Rb 

c  = 

— jsicC 

PB  = 

[  CBG  0  ]  0 

A  = 

F\q  4-  F2q  4-  F3A 

B.  Trimming  the  Aircraft 

Trimming  is  performed  for  both  zero  thrust  and  thrust  required  for  1-g  level  flight  based  upon  techniques 
outlined  in  Refs.  27  and  34.  A  cost  function  is  defined  as 


j  =  /Tf 


where  for  the  zero  thrust  or  gliding  cases 

f 


,  __  |  pitchi 


pitching  moment  about  the  origin  of  B  frame 
lift  —  weight 


(15) 


(16) 


For  the  case  of  1-g  level  flight,  the  longitudinal  B  reference  frame  linear  and  angular  accelerations  are  used, 
such  that 

f=  \  VB ,  }  (17) 
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The  cost  function  J  is  then  minimized  over  the  solution  space  using  the  elevator  deflection  angle,  <5e,  the 
body  angle  of  attack,  a,  and  thrust,  6t.  A  simple  numerical  Newton-Raphson  method  is  used  to  find  the 
local  minimum  of  the  search  variable,  i.e., 


(18) 


where 


(19) 


The  search  variable,  <S,  is  updated  by 


£fc+i  —  $k  4-  AS 


(20) 


and  fk+ 1  and  ^§s\k  are  recomPuted  using  5/ c+i.  The  process  continues  until  the  cost  function  J  reduces 

to  some  prescribed  tolerance.  To  prevent  divergence  of  the  solution,  Sk+ i  is  checked  at  each  iteration  step 
and  kept  within  a  prescribed  set  of  bounds.  The  Jacobian 


(21) 


is  computed  numerically  through  finite  differences.  The  entire  procedure  is  outlined  in  Figure  3. 

C.  Solution  of  EOM 

To  solve  the  nonlinear  differential  equations,  Eq.  1 4,  a  high  frequency  dissipative  time  stepping  approach  was 
implemented.  A  Modified  Newmark  Method  was  used  and  is  described  in  Ref.  35.  The  Modified  Newmark 
Method  was  selected  based  upon  its  ability  to  accurately  integrate  large  systems  of  equations  including  ones 
with  repeated  eigenvalues,  its  relative  ease  of  implementation  with  the  current  EOM  modeling,  and  the 
derivation  of  both  a  first  and  second  order  method.36,37 

III.  Trajectory  Control  of  Very  Flexible  Aircraft 

A.  Challenges  for  Trajectory  Control  of  Very  Flexible  Aircraft 

Due  to  the  nature  of  the  very  flexible  aircraft  construction,  the  wings  will  typically  have  a  lower  stiffness 
than  the  fuselage,  generating  stronger  coupling  of  the  rigid  body  and  structural  motion  in  the  lateral  axis 
than  in  the  longitudinal  one.  This  creates  a  variety  of  challenges  to  be  overcome  by  controllers.  The  first  is 
the  requirement  of  an  integrated  controller  which  handles  flexibility  as  well  as  rigid  body  motion.  Typically 
the  first  wing  bending  mode  of  this  class  of  aircraft  is  less  than  10  rad/s,  creating  direct  interaction  with 
classic  lateral  and  longitudinal  rigid  body  aircraft  modes  (spiral,  roll,  dutch  roll,  phugoid).  The  second  major 
challenge  is  a  time  delay  between  control  inputs  and  B  reference  frame  movement  due  to  the  flexibility  of 
the  aircraft.  This  time  delay  creates  a  non-minimum  phase  (NMP)  system  when  the  governing  differential 
equations  are  linearized.  The  third  major  challenge  is  the  introduction  of  additional  NMP  behavior  due  to 
adverse  yaw  from  aileron  inputs.  The  adverse  yaw  problem  may  be  avoided  through  the  use  of  spoilers  for 
roll  control,  but  that  is  not  addressed  in  this  work.  The  fourth  challenge  is  the  location  of  the  linearized 
structural  eigenvalues  near  the  imaginary  axis.  And  finally  the  generalized  mass  and  damping  matrices  are 
state  dependent. 

B.  Requirements  and  Assumptions  for  Trajectory  Control  of  Very  Flexible  Aircraft 

Before  designing  any  controller,  a  set  of  performance  objectives  should  be  established.  Currently,  despite 
operational  requirements,4  there  are  no  published  performance  specification  requirements  for  very  flexible 
aircraft  either  from  the  military,  civilian  government,  or  industry  authorities.  Given  this  void,  it  is  worthwhile 
to  refer  back  to  the  piloted  aircraft  military  standards,  MIL-STD-1797A  .38  While  very  flexible  aircraft  do 
not  have  the  safety  requirements  inherent  in  piloted  aircraft,  minimal  control  performance  requirements  are 
necessary  for  both  manned  and  unmanned  aircraft  to  complete  the  various  missions  and  potentially  fly  in 
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Figure  3.  Flow  for  Trim  Solution 
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the  National  Airspace  System.  For  classification  purposes,  a  very  flexible  vehicle  will  be  considered  a  large 
land  based  transport  type  aircraft,  Class  III-L,  as  in  Table  1.  From  MIL-STD-1797A  ,  Table  2  summarizes 


Table  1.  Aircraft  classes  as  defined  in  MIL-STD-1797A 


Class 

Type  of  Aicraft 

Example 

I 

Small  Light  Aircraft 

Primary  Trainer,  Light  Utility 

II 

Medium  Weight, 

Low-to-medium  maneuverability 

Search  and  Rescue,  I SR 

III 

Large,  Heavy, 

Low-to-medium  maneuverability 

Heavy  Transport,  ISR 

IV 

High  Maneuverability  Aircraft 

Fighter,  Attack 

Note:  a  -L  stands  for  land  based  and  a  -C  for  carrier  based 

the  desired  roll  requirements  for  Level  2  flying  qualities  (Table  3)  during  take-off,  climb,  loiter,  and  landing 
flight  conditions.  Furthermore,  the  maximum  bank  angle  required  shall  be  determined  by  the  bank  angle 


Table  2.  Final  roll  performance  objectives  for  very  flexible  aircraft 


Flight  Phase 

Bank  Angle  Change/ 
Seconds  to  achieve 

Roll  Mode  time  constant 

Climb  and  Loiter  -B* 

30°/3.9  s 

3.0  s 

Take-off  and  Landing  -C* 

3074.0  s 

3.0  s 

Note:  *  is  the  letter  designator  from  MIL-STD-1797A  for  the  particular  flight  phase 

Table  3.  Qualitative  degrees  of  suitability  and  levels  as  defined  in  MIL-STD-1797A 


Level 

Description 

1  -  Satisfactory 

Flying  Qualities  clearly  adequate  for 
the  mission  Flight  Phase 

2  -  Acceptable 

Flying  Qualities  adequate  to 
accomplish  the  mission  Flight  Phase 

3  -  Controllable 

Flying  Qualities  such  that  the 
aircraft  can  be  controlled  in  the 
context  of  the  mission  Flight  Phase  . . . 
workload  is  excessive  or  mission 
effectiveness  is  inadequate 

required  to  make  a  standard  two-minute  turn  (turn  rate  of  3°^  or  0.05236  ^).  Since  MIL-STD-1797A  does 
not  provide  specific  guidance  for  climb  rates,  a  maximum  climb  rate  of  2000  (10.16y)  at  sea-level  at 

maximum  gross  weight  shall  be  used.  This  rate  is  based  upon  reasonable  climb  rates  of  aircraft  similar  in 
size  and  weights. 

Three  basic  maneuvers  will  be  used  to  study  controller  performance.  The  first  maneuver  is  a  wings-level 
altitude  change.  The  second  maneuver  is  a  steady  level  turn  starting  from  a  zero  bank  angle.  And  the  third 
maneuver  is  climbing  turn. 

Finally  the  assumption  will  be  made  that  the  aircraft  will  be  flying  well  below  any  divergence,  flutter, 
or  limit  cycle  oscillation  boundaries.  While  this  assumption  is  restrictive  for  an  aircraft  dominated  by 
aeroelastic  effects,  it  is  an  important  first  step  in  the  development  of  controllers  for  performing  basic  aircraft 
maneuvering.  Future  research  should  focus  on  extending  the  controller(s)  beyond  aeroelastic  boundaries. 
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C.  Development  of  a  Control  Architecture  for  Very  Flexible  Aircraft 

As  discussed  in  Section  III. A,  very  flexible  aircraft,  and  particularly  HALE  ones,  present  unique  challenges 
to  the  design  of  a  controller.  This  section  details  the  specifics  of  a  proposed  very  flexible  aircraft  controller. 


1.  Traditional  Controller  Design  Difficulty 

Traditional  methods  of  modern  aircraft  control6,27,39  have  relied  upon  the  state  vector 
x  =  I  Vbx  Vbv  VBt  UJbx  Ubv  $0  4'  PBr  Pbu  Pbz 


(22) 


or  variations  of  x,  where  the  B  reference  frame  linear  and  angular  velocities  are  vb  and  ujb,  the  classic 
Euler  angles  of  roll  (positive  right  wing  down),  pitch  (positive  nose  up),  and  yaw  (positive  nose  right)  are 
given  by  $,  0,  and  $  respectively,  and  the  inertial  position  is  given  by  ps-  Linearizing  the  rigid  body 
governing  differential  equations  about  this  state  vector,  and  determining  a  constant  gain  matrix,  A',  has 
been  shown  in  numerous  papers  and  books  to  yield  satisfactory  results  when  applied  to  nonlinear  flight 
dynamic  models.  However,  the  current  problem  with  additional  aeroelastic  effects  has  rendered  this  to  be 
ineffective  for  trajectory  control.  In  the  process  of  developing  a  stabilizing  controller  architecture  for  very 
flexible  aircraft,  traditional  controller  designs  were  initially  applied  to  a  statically  deformed  rigid  aircraft 
model.  This  controler  architecture  was  shown  to  have  difficulties  in  providing  closed  loop  stable  trajectory 
tracking.33 


2.  A  Heuristic  Approach  Mimicking  a  Human  Pilot 

Due  to  the  difficulty  in  finding  a  stable  controller  for  the  rigid  body  using  traditional  techniques,  the  method 
was  not  attempted  with  additional  elastic  states.  A  method  of  decoupling  the  linear  and  nonlinear  effects 
of  the  aircraft  response  was  designed  based  loosely  upon  aircraft  pilot  training.  A  well  trained  human  pilot 
has  been  taught  to  command  flight  path  angle,  7,  bank  angle,  /i,  and  their  rates,  7  and  //,  when  changing 
altitude  or  heading  4/,  through  a  fast-  (inner)  and  slow-  (outer)  loop  process.  The  fast- loop  consists  of 
commanding  pitch  rate,  pitch  angle,  roll  rate,  and  sideslip  angle.  These  changes  are  commanded  by  the  pilot 
typically  through  the  use  of  four  controls:  elevator  (£e),  aileron  ( Sa ),  rudder  (6r),  and  throttle  (<5t).  Once 
an  angle  and  its  rate  have  been  satisfactory  set,  the  pilot  then  perforins  a  slow-loop  (outer-loop)  function  of 
cross  checking  altitude,  rate  of  climb,  and  other  states.  Heffley  et.  al.40  have  described  this  type  of  inner- 
and  outer-loop  control  as  an  implicit  inner-loop  and  explicit  outer-loop  tracking  task.  Additionally  for  level 
turns,  altitude  changes,  and  landing  flair  maneuvers,  Heffley  and  Schulman41  have  provided  pilot  models  for 
outer-loop  control. 

In  similar  fashion,  this  work  proposes  building  a  stable  inner-loop,  utilizing  both  linear  and  nonlinear 
theory,  where  the  states  are  the  linear  and  angular  B  reference  frame  velocities,  /i,  Eq.  5,  and  augmented  error 
states.  The  outer-loop  or  slow-loop  is  then  controlled  through  a  nonlinear  transformation  and  traditional 
proportional,  integral,  derivative  (PID)  control  techniques.42  Conceptually  the  entire  system  is  shown  in 
Figure  4.  In  the  figure  /Jo  are  the  initial  trimmed  B  reference  frame  linear  and  angular  velocities.  For  the 
outer-loop  a  commanded  B  reference  frame  velocity,  /JCOm  and  differences  are  generated 

A[3Com  =  ficom  ~  ft  (23) 

and  used  to  excite  inner-loop  integrator  error  states.  For  the  inner-loop  the  difference  in  (3 


A  U  =  U-Uo  (24) 

is  required  to  generate  differences  in  the  control  vector  from  the  initial  control  deflections,  uo,  such  that 

u  =  u o  4-  A u  (25) 

and  An  is  a  function  of  the  gain  matrix,  A',  and  the  augmented  error  states,  e 

Au  =  K(A/3,e)  (26) 

Also  the  subscripts  des  and  com  are  the  desired  and  commanded  values. 
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Outer-Loop  (Slow-Loop) 


Figure  4.  Controller  Concept  Mimicking  a  Human  Pilot 


3.  Nonlinear  Transformation  from  Flight  Path  and  Bank  Angles  to  Body  Velocities 

The  classic  rigid  body  aircraft  kinematic  equations  used  for  the  nonlinear  transformation  are2  '■ 39,43 


\JlL2  +V2  +  W2 

(27) 

s\n~l(v/V) 

(28) 

sin-1  ( (cos  a  sin  0  -  sin  a  cos  0  cos  4>)  cos/3s  —  cos0sin4>sin^s) 

(29) 

,  /  cos  0  sin  $  cos  Ba  4-  (cos  a  sin  0  —  sin  a  cos  0  cos  4>)  sin  (3a\ 
sm  M 

\  cos  7  / 

(30) 

P  4  tan  0  ( Q  sin  <1 >  +  P  cos  4>) 

(31) 

Q  cos  4>  -  R  sin  4> 

(32) 

where  V  is  the  total  airspeed  and  the  classic  aircraft  longitudinal,  lateral,  and  vertical  velocity  components 
are  u  positive  out  the  nose,  v  positive  out  the  right  wing,  and  w  positive  down.  Further  7  is  the  flight  path 
elevation  angle,  a  is  the  aircraft  angle  of  attack,  p  the  flight  path  bank  angle,  and  fja  the  sideslip  angle. 
The  classic  aircraft  angular  rates  are  roll  rate,  P,  pitch  rate,  <2,  and  yaw  rate,  R.  The  dynamic  equations  of 
interest  are 

L  =  ±pV2SCLaa  (33) 

/  F 

v  =  —  Ru  4-  Pw  4-  gn  sin  <I>  cos  0  H — -  (34) 

m 

where  L  is  the  aircraft  lift,  p  is  the  atmospheric  density,  5  is  the  surface  area  of  the  wing,  Cl0  is  the 
equivalent  aircraft  lift  curve  slope,  g0  is  the  magnitude  of  the  gravity  vector,  Fy  is  the  lateral  force  due  to 
aerodynamic  and  control  inputs,  and  m  is  the  aircraft  mass.  Equation  33  is  a  steady  state  lift  approximation 
and  will  be  used  to  develop  an  angle-of-attack,  a,  dynamic  relationship.  Equation  34  is  the  rigid  body  lateral 
acceleration  EOM  and  will  be  used  to  develop  relationships  for  the  B  reference  frame  angular  velocities.  The 
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angle  of  attack  and  sideslip  are  related  to  the  longitudinal  velocites  as 


u 

f  V'cos(a)  COs(jtfs) 

<  V 

►  =  <  Vsin(A) 

w 

1  1/ sin(a)  cos(/*g) 

such  that  the  angle  of  attack  may  be  written  as 

a  =  tan-1  (w/u) 


(35) 


(36) 


Before  Eqs.  27  -  86  can  be  used  in  the  current  formulation,  transformations  of  the  linear  and  angular  body 
velocities  and  quaternion  to  classic  Euler  angles  are  presented.  For  the  linear  and  angular  velocities  the 
relationships  are 


and 


f  > 

VBX 

"  0 

1 

0 

t 

u 

'  vBy 

>  = 

1 

0 

0 

< 

v  > 

Vb x 

_  0 

0 

-1 

w 

UBX 

’01  0 

p  ' 

Wfl, 

►  = 

1  0  0 

Q 

0 

0 

1 

1 _ 

R 

\ 

(37) 


(38) 


To  develop  the  quaternion  relationship,  first  the  body,  B ,  to  inertial,  (7,  rotation  matrix  in  terms  of  the 
Euler  angles  is  found  as 


(c4>  c'P  4-  s4>  s0  5^) 

(— cty  4-  sty  s0  cty) 

—  sty  C0 

C0  s'!/ 

C0  cty 

S0 

(sty  cty  -  cty  S0  sty) 

—  (sty  sty  4  c4>  s0  c^) 

cty  C0 

(39) 


where  the  shorthand  costy  =  cty  and  sinty  =  sty  is  used.  The  Euler  angles  of  interest  have  the  following 
relationship  to  the  quaternions33 


sinB 

=  2(C2C3  +  <oCl) 

(40) 

cos  0  cos  ty 

=  Co  -  C?  -  <2  +  Cl 

(41) 

cos  0  sin  ty 

=  2(— C1C3  +  C0C2) 

(42) 

To  perform  the  transformation  from  given  inputs  of  flight  path  angle  and  bank  angle  to  desired  body 
velocities,  the  following  assumptions  are  made: 

1.  Angle  of  side-slip,  (j8,  is  zero; 

2.  The  Euler  bank  angle,  ty  will  be  used  in  place  of  bank  angle  //; 

3.  Total  velocity,  V,  is  prescribed  and  typically  assumed  constant; 

4.  Angle  of  attack,  a,  is  proportional  to  costy; 

5.  Lateral  side  forces  are  only  a  function  of  cross  coupling  of  linear  and  angular  velocities  and  gravity 
component  due  to  Euler  bank  and  pitch  angle,  and  lateral  velocity  is  constant. 

While  zero  side  slip,  assumption  1,  is  typically  a  requirement  of  manned  aircraft  for  comfort  reasons,  it  is 
also  desired  to  minimize  coupling  of  the  lateral  and  longitudinal  aircraft  motion.  Assumption  2  is  justified  by 
applying  small  angle  assumptions  to  7,  Euler  pitch  angle,  0,  and  angle  of  attack,  a.  Using  these  assumptions 
and  Eq.  2*) 

H  «  ty  (43) 

Further,  if  a  navigational  loop  is  wrapped  around  the  architecture  described  in  Figure  1,  the  Euler  roll  angle 
would  be  convenient  as  an  input.  Assumption  3  is  made  to  simplify  the  resulting  nonlinear  transformation 
equations.  Assumption  4  has  two  significant  effects.  The  first  is  unsteady  aerodynamic  effects  presented 
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in  Ref.  32,  are  neglected.  This  is  done  because  one  of  the  goals  of  this  work  is  to  development  an  initial 
control  architecture  for  very  flexible  aircraft,  without  consideration  of  aeroelastic  boundaries.  Additionally 
the  controller  is  designed  to  maintain  aircraft  velocities  well  below  these  aeroelastic  boundaries,  preventing 
destabilization  of  the  controller  through  unmodeled  dynamics  (within  the  controller).  The  second  significant 
effect  of  assumption  4  is  an  increase  in  angle-of-attack  to  generate  additional  lift  for  a  wings-level  climb  is  not 
considered.  This  is  reasonable  approximation  utilized  in  most  rigid  body  aircraft  performance  equations43 
where 

D  .  r^r  u  Excess  Thrust 

Rate  of  Climb  = - .  - -  (44) 

Weight 

Assumption  4  is  derived  from  the  steady  level-turn  performance  equation43 


Tlload  — 


1 


cos4> 


(45) 


where  the  load  factor,  n/oad,  is  the  non-dimensionalized  acceleration  of  gravity.  Assuming  linear  three 
dimensional  lift  theory  and  a  level  turn 

L<t>  =  ni0adL^= o  (46) 

Using  Eq.  33  in  Eq.  16,  it  can  be  shown  that 

a#  = - —ao  (47) 

cos4> 

where  Qo  is  the  angle  of  attack  corresponding  to  a  wings  level  at  steady  level  flight.  Finally  assumption  5, 
which  simplifies  the  lateral  EOM,  Eq.  34,  provides  one  of  three  resulting  equations  used  to  determine  the  B 
reference  frame  angular  velocities. 

Using  Eqs.  27  -  47  and  the  five  assumptions,  the  resulting  set  of  equations  are: 


•  angle  of  attack  relationships: 

a  = 

1 

cos<f> 

(48) 

a  = 

tan  4>  sec  4>  ao  $ 

(49) 

•  linear  body  velocity  relationships: 

VBZ 

=  —V  sin  q 

(50) 

VBZ 

=  —  V  cos  a  a 

(51) 

vBy 

=  v^2  -  4 

(52) 

vbv 

VBZVBZ 

VB„ 

(53) 

•  flight  path  angle,  Euler  pitch  angle,  and  angle  of  attack  relationships: 


sin  7 
0 


cos  a  sin  0  —  cos  $  sin  a  cos  0 

cos  77  4-  (sin  a  sin  0  +  cos  4>  cos  a  cos  0)  d  ^ 

cos  a  cos  0  4-  cos  4>  sin  a  sin  0 
sin  4>  sin  q  cos  0  <t> 
cos  a  cos  0  4  cos  4>  sin  a  sin  0 


(54) 


(55) 


•  body  angular  velocity  relationships. 


4> 


g0  sin  4>  cos  0 
0 


(tan  0  sin  4>)  ujbx  4  ujbv  —  (tan  0  cos  4>)  ujbx 
vb 2  ubv  ~  vbu  ub2 
cos  <t>  ujbs  4  sin  4>  ubz 


(56) 

(57) 

(58) 
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Given  the  flight  path  angle  and  rate,  7  and  7,  and  the  Euler  roll  angle  and  rate,  <I>  and  4>,  Eqs.  IS  -  58  are 
used  to  solve  for  the  commaned  B  reference  velocities,  f3com. 

To  close  the  outer-loop  of  Figure  I,  the  actual  flight  path  angle,  7,  Euler  bank  angle,  3>,  and  their  rates, 


7  and  <I>  are  determined. 

Differences  are  then  found 

A 7  =  7-  7 des 

A7  =  7  -  7 des 

A<t>  =  <t>  -  $des 

A<t>  =  <i>-  <j><ies 

(59) 

The  commanded  flight  path  elevation  angle  and  Euler  roll  angle  are  computed  using  a  traditional  linear 
technique  of  proportional,  integral,  derivative  (PID)  control 

7com  =  7 des  "1"  A  7 com 

ifcom  =  'Idea  4"  Alcorn 
^ com  =  ^ des  “1“  A *&com 

^com  =  ^ des  4-  A^coro 

(60) 

where 

A7  com 

=  k~,Ay  +  kiy  J  (A-y)  dr  +  kny  J  (Ay)  dr)  dr  +  k^  Ay 

A' 7com 

=  k^Ay  +  fc/->  A7  +  kui  ( A7)  dr  +  fc///-y  J  (Ay) dr'j  dr 

^^com 

=  4A4>  4-  ki$  f  (A<I>)  dr  4-  kd* A3> 

Jo 

(61) 

^^com 

=  **Ai  +  k7*A*  +  fc/7*  [  (A4>)  dr 

The  outer-loop  feedback  gains,  fc7,  fc/7,  etc,  of  Eq.  61  are  initially  determined  using  the  SISO  heuristic 
guidelines  of  Ziegler  and  Nichols44,45  where 


u  =  kce  4-  ki  /  e  dt  4-  kje 


/■ 


(62) 


For  PID  and  PI  controllers,  Ziegler  and  Nichols  recommend  values  are  given  in  Table  1,  where  ku  is  the 
closed  loop  gain  required  to  make  the  system  marginally  stable  and  tu  is  the  corresponding  period  between 
oscillations.  The  gains  are  then  adjusted  to  meet  desired  performance.  Due  to  the  potential  instability  of 


Table  4.  Ziegler  and  Nichols44, 46  PID  and  PI  tuning  parameters 


PID 

PI 

Proportional  Gain 

kc  =  0.6/cu 

kg  —  0.45fcu 

Integral  Gain 

=  05tT 

h  r  —  _ fin - 

A,/  0.85t„ 

Derivative  Gain 

kd  — 

integral  feedback,42  the  additional  integral  gains  (fc//7,  fc//7,  fc///7,  etc.)  are  chosen  to  be  at  least  an  order 
or  more  of  magnitude  less  than  the  first  integral  gains  ( kj1 ,  fc/7,  etc.). 
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D.  Modification  of  Proposed  Control  Architecture  to  a  Flexible  Vehicle 

This  section  presents  a  method  for  extending  the  control  architecture  previously  introduced  to  a  flexible 
aircraft.  First  a  separation  of  the  lateral  and  longitudinal  motion  is  presented.  Then  control  techniques 
are  applied  to  the  lateral  and  longitudinal  motion  separately.  Finally  coupling  between  the  lateral  and 
longitudinal  motion  is  shown  to  be  handled  in  an  outer-loop  strategy. 

Section  III.C  outlined  a  control  strategy  for  very  flexible  aircraft.  The  philosophy  of  this  approach  is  to 
provide  a  systematic  method  where  the  dominant  coupling  of  lateral  and  longitudinal  motion  is  handled  in  a 
slower  outer-loop,  and  B  reference  frame  linear  and  angular  velocities  are  handled  in  a  faster  inner-loop.  The 
six  states  (linear  and  angular  B  reference  frame  velocities)  of  the  inner-loop  are  typically  easily  separated  into 
longitudinal  and  lateral  motion  with  corresponding  control  effectors.  This  allows  individual  control  schemes 
to  be  applied  to  each  set  of  the  inner-loop  dynamics  (longitudinal  and  lateral).  The  outer  loop  controls  the 
required  longitudinal  and  lateral  motion  necessary  for  trajectories  where  this  motion  is  coupled,  i.e.  steady 
level  turns  and  climbing  turns.  Additionally  the  outer-loop  handles  coupling  of  the  lateral  and  longitudinal 
motion  due  to  gravitational  effects.  While  this  idea  is  not  new  and  has  been  applied  to  high-angle-of-attack 
flight,46'48  it  typically  has  not  been  applied  to  “low-  to  moderate-angle-of-attack  flight  regimes”.47  This  is 
due  largely  to  traditional  aircraft  control  design  assuming  relatively  small  to  moderate  amounts  of  nonlinear 
cross  coupling.  With  very  flexible  aircraft,  the  large  potential  movement  of  the  cm  from  the  origin  of  the  B 
reference  frame  can  create  a  significant  nonlinear  coupling. 

1.  Separation  of  Lateral  and  Longitudinal  Motion 

As  shown  by  Shearer  and  Cesnik, 32,33  for  this  class  of  vehicles  with  relatively  stiff  fuselages,  longitudinal 
motion  does  not  appear  to  be  significantly  affected  by  wing  flexibility.  So  it  is  assumed  that  wing  flexibility 
is  a  secondary  and  minimal  contribution  to  aircraft  longitudinal  motion.  Using  this  assumption  and  the 
previous  assumptions  of  Section  111.0.3,  a  separation  between  longitudinal  and  lateral/elastic  motion  is 
made.  Additionally  an  examination  of  the  rigid  body  LTI  state  matrices,  Ref  33,  supports  this  assumption 
of  minimal  cross  coupling  between  the  lateral  ( vbt ,  ujbv,  and  ljbz)  and  the  longitudinal  ( vbu ,  vst,  and  u ;#x) 
states.  An  eigenvalue  analysis  of  a  representative  HALE  aircraft  further  supports  this  assumption. 

The  lateral  and  the  longitudinal  eigenvalues  and  eigenvectors  of  the  B  reference  frame  for  a  statically 
deformed  aircraft  at  a  given  (heavy)  fuel  weight  are  presented  in  Tables  5  and  0.  As  it  can  be  seen  by 
examining  the  eigenvectors,  there  exists  a  fairly  distinct  separation  of  lateral  and  longitudinal  motion. 
When  the  elastic  states  are  included,  a  distinct  change  in  the  lateral  eigenvalues  and  eigenvectors  occurs, 
as  seen  in  Table  7.  The  “Sideslip”  and  “Sideslip/Yaw  Rate”  modes  experience  a  significant  change  in  the 
eigenvalues  and  changes  in  the  sign  of  the  eigenvector  components.  However  the  real  change  is  in  the  “Roll” 
mode  as  it  no  longer  has  any  significant  contribution  from  the  rigid  body  roll  rate,  u/b„.  But  rather  u jbu 
now  contributes  to  various  elastic  strains  and  rates  (not  presented).  Because  of  this,  the  elastic  states  are 
assumed  to  be  tightly  coupled  with  the  rigid  body  lateral  motion.  The  longitudinal  eigenvalues  and  vectors 
also  experience  a  change  as  seen  in  Table  s,  with  a  major  change  in  the  “Vertical/Longitudinal”  eigenvalue. 
This  eigenvalue  now  has  contributions  of  about  2%  from  several  of  the  in  and  out  of  plane  bending  strain 
rates  (not  presented).  However  from  Table  8,  over  94.5%  of  the  contribution  comes  from  the  B  reference 
frame  longitudinal  states.  Because  of  this  significant  contribution,  it  assumed  that  the  longitudinal  states 
are  decoupled  from  the  elastic  states.  While  this  strong  coupling  of  lateral  and  elastic  modes  cannot  be 
generalized  for  all  very  flexible  aircraft,  the  analysis  does  hold  for  the  representative  aircraft  used  here  at 
different  loadings  (fuel)  conditions  (not  shown).  In  general,  an  eigenvalue  and  eigenvector  analysis  should 
be  performed  at  nonlinear  equilibrium  conditions  of  a  very  flexible  aircraft  configuration  to  determine  the 
coupling  of  elastic  states  with  B  reference  frame  motion.  Based  upon  the  outcome  of  the  analysis,  various 
linear  and  nonlinear  control  schemes  could  be  utilized.  For  example,  Gregory22"25  has  developed  filtered 
nonlinear  control  techniques  and  optimum  sensor  placement  for  longitudinal  control  with  coupled  aeroelastic 
effects.  In  this  dissertation  a  linear  LQR  formulation  is  developed  for  the  lateral  motion  and  a  nonlinear 
dynamic  inversion  for  the  longitudinal  motion. 

2.  Lateral  Motion  Inner-Loop  Controller 

Due  to  the  inherent  NMP  zeros  present  in  the  lateral  dynamics,  Section  II LA,  most  nonlinear  control  schemes 
are  not  sufficient  for  controlling  it.  Because  of  this,  a  traditional  LQR  controller  is  used  for  the  inner-loop 
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Table  5.  Lateral  eigenvalues  and  eigenvectors  for  rigid  body  aircraft  with  elastically  deformed  wing 


Name 

Eigenvalue  (rad/s) 

Normalized  Eigenvector 

Sideslip 

-2.5081  •  10~4 

9.9729  •  10"1 
-4.0745  •  10“12 
5.4152  •  10"13 
2.9881  •  10-14 
-2.3031  •  10~3 
-7.3474  •  10~2 

VBX 

VBv 

VBX 

Vbx 

VBv 

VBr 

Sideslip/ Yaw  Rate 

—  1.0132 

9.9003  •  10"1 
1.9776  •  10“15 
5.4550  •  10"15 
1.8435  •  10“14 
-5.8755  •  10“2 
1.2801  •  10 1 

VBX 

VBv 

VBZ 

Vbx 

vbu 

UB* 

Roll  Rate 

-6.7172 

2.9039  •  10“2 
-3.4145  •  10~13 
9.5019  •  10“12 
-8.7481  •  10-13 
9.9957  •  10-1 
4.8056  •  lO"3 

VBX 

VBv 

VBZ 

VBX 

V By 

Vbz 

Table  6.  Longitudinal  eigenvalues  and  eigenvectors  for  rigid  body  aircraft  with  elastically  deformed  wing 


Name 

Eigenvalue  (rad/s) 

Normalized  Eigenvector 

2.2795 

IQ-14 

VBX 

1.1025 

io-2 

VBy 

Vertical 

-3.1402 

-9.9077 

10-1 

VBZ 

1.3514 

10-* 

VBX 

5.5100 

10-13 

VBv 

6.9795 

10-15 

UJbz 

5.3907 

10"13 

VBX 

-5.9205 

10-1 

VBy 

Vertical/Longitudinal 

-7.0227  •  10“2 

8.0334 

10"1 

^Bz 

6.4225 

10  2 

VBx 

-1.9984 

10-13 

VBy 

-3.3784 

10~14 

VBZ 

2.2961 

10-13 

VBX 

9.7330 

10-1 

VBv 

Longitudinal 

1.8195  •  10"2 

2.2882 

10-1 

VBZ 

1.8057 

10  2 

VBX 

2.3820 

10-14 

VBy 

-1.7343 

10"14 

vbz 
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Table  7.  Lateral  eigenvalues  and  eigenvectors  for  flexible  aircraft 


Name 

Eigenvalue  (rad/s) 

Normalized  Eigenvector 

Sideslip 

-1.5489  •  10-1 

-9.9566  •  10"1 
-6.4196  •  10"13 
-1.8588  •  10-13 
-2.3153  •  10“14 
4.7108  •  10“2 
-8.0219  •  10"2 

VBX 

VBv 

VBZ 

UBX 

V By 

UBt 

Sideslip/ Yaw  Rate 

5.3050  •  10"2 

-9.9163  •  lO-1 
7.0617  •  10“12 
1.7969  ■  10"12 
9.3405  •  10"14 
-1.8629  •  10~2 
1.2773  •  lO-1 

VBX 

VBy 

vBt 

^Bx 

UBV 

ubz 

Table  8.  Longitudinal  eigenvalues  and  eigenvectors  for  flexible  aircraft 


Name 

Eigenvalue  (rad/s) 

Normalized  Eigenvector 

-1.2290  • 

10-12 

VBX 

2.0530  • 

10-1 

VBv 

Vertical 

-2.5686 

-9.7172  • 

10"1 

VBZ 

7.4426  • 

10-2 

^Bx 

-3.2549  • 

10-13 

UBV 

-3.8281  • 

10-13 

ubz 

3.0616  ■ 

10-12 

vbx 

-2.5808  • 

10"1 

VBv 

Vertical/Longitudinal 

-3.2479 

9.0910  • 

10'1 

VBZ 

2.3189  • 

10  2 

ubx 

1.2535  • 

10-12 

UBV 

9.3817  • 

10-13 

UBZ 

-5.5698  • 

10-12 

VB, 

-9.7008  • 

10-1 

VBv 

Longitudinal 

4.7465  •  10-2 

-2.4244  • 

10_1 

VB, 

-1.2841  • 

10-2 

WB, 

-8.6632  • 

10-14 

UBy 

6.6309  • 

10-13 

Wfl, 

18  of  IX 


American  Institute  of  Aeronautics  and  Astronautics 


lateral  motion.  For  the  lateral  controller,  full  state  feedback  was  assumed,  where  the  states  to  be  controlled 
are  ^ 

Xiat  =  [  eT  eT  vBt  vBy  uBt  ]  (63) 

While  the  assumption  of  full  elastic  state  feedback  (e  and  e)  may  present  a  practical  limitation,  it  is  used 
here  as  a  starting  point  for  very  flexible  aircraft  control  architecture  development.  Linear  time  invariant  A 
and  B  matrices  are  generated  from  Eq.  I  l  for  the  lateral  states  of  Eq.  (>3.  Error  states  are  then  augmented 
to  the  system  such  that 

Xiat  =  [  eT  eT  vBx  U )By  UE 

and 

vBx  -0 

^ By  —  L^ByCOm 

k  UJBt  -U)BzCOm 

The  augmented  system  can  then  be  represented  as 


e/at  =  < 


lat  — 


Alat  0 
Clat  0 


xlat  "I" 


Blat 

0 


Ulat  + 


0 

-I 


0 

^  By  com 
k  ^ Bzcom 


where 

’  0  0  1  0  0 

Clat  =  0  0  0  1  0 

0  0  0  0  1 

The  lateral  control  vector  is  found  using  the  control  law 


Ulat 


= -  [  } 


(64) 

(65) 


(66) 


(67) 


(68) 


where  Kiat  and  Keiat  are  found  using  standard  LQR  techniques  applied  to  the  augmented  system  of  Eq.  <i0. 


3.  Longitudinal  Motion  Inner-Loop  Controller 

In  a  similar  manner  to  Al-Hiddabi49  and  Al-Hiddabi  and  McClamroch,50  partial  feedback  linearization  or 
dynamic  inversion  is  used  for  development  of  the  very  flexible  aircraft.  Here  dynamic  inversion  is  applied  to 
the  longitudinal  states  of  interest,  i.e., 


Xiong  =  [  t >BU  VB,  UBr  ] 

Dynamic  inversion51,52  is  a  process  by  which  nonlinear  differential  equations  of  the  form 

x  =  f(x)  +  g(x)u 
V  =  h{x) 

are  transformed  into  an  equivalent  linear  form 

X=v 

The  transformation  is  accomplished  by  first  making  a  change  of  variables 

x  =  y 

and  then  finding  \ 

.  .  dh  . 

x=y  = 

dh  dh 

=  tef{x)  +  te9{x)u 


(69) 

(70) 

(71) 

(72) 


(73) 
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A  new  control  variable,  v,  is  then  defined  as 


dh  dh 
v  =  d^f(x)  +  d^g{x)u 


(74) 


and  finally  solving  for  u 


dh  ,  \ 
^9{X) 


(75) 


The  transformed  system,  Eq.  7 1 ,  can  then  be  augmented  with  error  states  if  desired,  and  the  control  variable, 
v  is  solved  as 


V  =  -K\ 


(76) 


using  any  number  of  linear  controller  techniques,  yielding  closed  loop  trajectory  tracking. 

For  the  very  flexible  aircraft,  a  subset  of  the  linear  and  angular  body  velocities,  /?,  are  desired  to  be 
controlled,  Eq.  (i!»,  using  dynamic  inversion.  To  accomplish  this,  0  is  found  from  the  EOM  for  e  and  (3 . 
Starting  with  Eq.  I  I, 


Alp'f'  Mfb 
Mbf  Mbb 


Cff  Cfb  f  *  1  _  f  Kff e  1  +  /  Rf  1 
Cbf  CBb  \  P  J  \  0  J  |  ^  j 


(77) 


these  can  be  solved  by  inverting  the  generalized  mass  matrix.  Using  Fact  2.15.3  of  Ref.  53 

A  B 
Bt  C 

(A-  BC~lBT)~l  -(A-  BC~lBTy1  BC~l 

—C~1BT  ( A  -  BC~iBT)-'  C~lBT  (A  -  BC~lBT)~1  BC~l  +  C~l 

and  noting  Mfb  =  Mbf  from  Eq.  7,  the  inverse  of  the  generalized  mass  matrix  of  Eq.  77  is 


Mff  Mfb 
Mbf  Mbb 


where 

Qm  =  {Mff  ~  MfbMqqM Bf)  (SO) 

Therefore  e  and  0  are  given  by 


e  =  -C\\i  -  C\20  -  QmKff*  4  Rfi 

(81) 

0  =  —  C220  —  C21I  4  Rbi 

(82) 

where 

C\\ 

=  QmCff  ~  QmMFbMbbCbf 

C\2 

=  QmCfB  -  Q M M F bMbbC SB 

C21 

=  ~MbbMBfQmCff  4  ( MbbMbfQmMfbMbb  4  MBB)  Cbf 

(83) 

C22 

=  - MBBM B fQ mC f b  4  ( MbbMbfQmMfbMbb  4  MBB)  Cbb 

Rf  1 

=  Qm  Rf  ~  QmMFbMbbRb 

Rbi 

=  ~MbbMBfQmRf  4  {MbbMbfQmMFbMbb  4  MBB)  RB 

The  generalized  force  vectors,  Rf  and  Rb,  can  be  further  expanded  into  a  control  affine  form 


Qm 

-MqXbMbfQm 


-Qm  MFbM  BlB 

MbbMBfQmMFbMbb  +  Mbb 


(79) 


Rf  =  Rfu= o  (e,  0,  C  A)  4-  RFu  (c,  €,  0,  <,  A )u  (84) 

Rb  =  Rbu=0(6,€,^,C,A)  4 /?Bti(£,6,^,C,A)u  (85) 
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The  body  velocity,  Eq.  *2,  can  then  be  written  in  the  generalized  form 

P  =  fp  +  gpu  (86) 

where 

f(3  =  -C22P  ~  C2\e  +  -MbbMbfQmRfu=0  + 

+  (MbbMbfQmMfbMbb  +  M bb )  Rbu=0  (87) 


and 

9(3  —  -Mb1bMbfQm Rfu  +  {MbbMbfQmMfbMbb  +  M BB )  Rbu  (88) 

Stability  of  dynamic  inversion  requires  the  so-called  internal  or  zero  dynamics54  to  be  stable.  For  the  partial 
Dynamic  Inversion  scheme  proposed,  the  internal  dynamics  consist  of  the  controlled  lateral  states  and  the 
uncontrolled  inflow  states.32  The  assumption  is  made  that  the  internal  dynamic  states  are  stable.  This 
assumption  is  valid  provided  the  linear  controller  of  Section  1 I  I.D.2  yields  closed  loop  stability  and  the  inflow 
states,  A,  are  not  driven  unstable  due  to  the  aeroelastic  boundaries  of  flutter  and  divergence.  While  avoiding 
the  aeroelastic  boundaries  is  a  limitation  of  the  current  design,  it  is  a  necessary  requirement  in  the  developing 
stages  of  controller  design  for  very  flexible  aircraft. 

The  longitudinal  outputs  to  be  tracked,  Eq.  6fl,  can  be  written  as 


where 

The  linearized  form  is 

where 


Vlong  —  C long  P 


Clong  ~ 


0  1  0  0  0  0 
0  0  1  0  0  0 
0  0  0  1  0  0 


Xiong  —  Vlong 
Vlong  —  9  p  {.Vlong  //?) 


The  system  is  then  augmented  with  longitudinal  error  states 


(89) 


(90) 


(91) 

(92) 


{Vbv  VBycom 
VBt  -  VBzcom  > 
UBX  ~  UBrcom  , 


and  the  final  linearized  system  is 


Xiong  1 

I 

{Vlong}  + 

0 

< 

*  ' 

V  ByCom 

VBxcom  * 

Clong  J 

0 

-/ 

K  & BTcom  4 

The  control  variable,  viong,  is  then  computed  using  full  state  feedback  as 


Vlong 


~[ 


Ky 


(93) 


(94) 


(95) 


where  the  gains,  KXlong  and  Keiong ,  are  computed  using  the  LQR  technique.  The  actual  control  signal,  uiong 
is  computed  using  Eq.  02. 

Finally,  since  dynamic  inversion  becomes  unstable  with  the  introduction  of  NMP  zeros,55  flight  path  angle 
is  tracked,  instead  of  altitude,  in  the  outer-loop  design  to  prevent  the  introduction  of  NMP  zeros  into  the 
inner-loop  longitudinal  dynamics.  This  is  necessary  since  in  traditional  wing/body/tail  aircraft  an  increase 
in  altitude  is  accomplished  by  deflecting  the  trailing  edge  surface  of  the  elevator  upwards.  This  creates  a 
downward  force  which  temporarily  decreases  the  overall  aircraft  lift  and  a  subsequent  decrease  in  altitude. 
However,  this  downward  force  also  creates  a  nose  up  pitching  moment,  which  increases  the  angle  of  attack. 
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Once  the  airflow  adjusts  to  this  increase  in  angle  of  attack,  overall  lift  is  increased  and  the  aircraft  begins 
to  climb.  In  mathematical  terms,  this  behavior  represents  a  NMP  system.  To  control  altitude  a  first-order 
approximation  of  the  rate  of  climb,  h,  is  used 


h  =  V  7  (96) 

By  prescribing  the  altitude  and  total  velocity,  V,  trajectories,  h,  7,  and  7  can  be  commanded  to  the  outer- 
loop. 


IV.  Trajectory  Control  Studies 

In  this  section,  trajectory  control  studies  are  performed  on  the  proposed  controller  architecture.  The 
control  architecture  is  initially  applied  to  an  elastic  aircraft  at  a  heavy  weight  condition.  Simulated  climbs, 
bank  angle  changes,  and  climb  and  turning  trajectories  are  presented  and  discussed.  The  controller  is  then 
applied,  without  changes  to  any  of  the  gains,  to  a  light  weight  condition  which  has  a  commanded  climb  and 
turn  trajectory.  Finally  a  representative  mission  profile  segment  is  developed  and  simulated. 

A.  Representative  HALE  Aircraft 

For  the  aircraft  simulation  studies,  a  representative  HALE  type  aircraft  was  created  and  it  is  shown  in 
Figure  5.  The  relevant  physical  properties  are  summarized  in  Table  9.  The  vehicle  was  designed  to  be 


Figure  5.  Representative  HALE  aircraft  model  (units  are  in  meters) 

statically  stable  for  moderate  wing  deflections  in  both  the  longitudinal  and  lateral  axes.  Table  10  also 
includes  the  trimmed  longitudinal  controls  and  state  (elevator  angle,  thrust  level,  angle  of  attack)  for  both 
a  gliding  and  1-g  level  flights.  Additional  details  of  the  vehicle’s  mass,  structural  damping,  and  stiffness 
parameters  are  provided  in  Ref.  33. 

The  vehicle  is  a  conventional  wing/body/tail  configuration  with  twin  vertical  tails.  It  is  representative  of 
a  HALE  aircraft  concept  being  considered  by  the  USAF.  The  aircraft  has  the  conventional  control  surfaces 
of  elevator,  aileron,  and  twin  rudders.  The  elevator  is  such  that  a  positive  elevator  control  input,  de,  results 
in  a  negative  pitching  moment  (nose  down).  The  left  and  right  ailerons  have  a  -1:1  gearing  ratio,  such  that 
a  positive  aileron  control  input,  <)'a,  results  in  a  roll  to  the  left  (left  wing  down).  The  twin  rudders  have 
a  1:1  gear  ratio  such  that  a  positive  rudder  control  input,  Jr,  produces  a  positive  yawing  moment  (nose 
left).  Recall  the  B  reference  frame  orientation  is  x  positive  out  the  right  wing,  y  positive  out  the  nose,  and 
z  positive  up.  Thrust  is  accomplished  using  a  simple  point  force  applied  at  the  origin  of  the  B  reference 
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Table  9.  Geometric  properties  of  the  very  flexible  aircraft  model 


Model  Parameters 

Property 

Value 

Units 

Fuselage  Length 

26.4 

m 

Wing  Span 

58.6 

m 

Wing  Area 

196.3 

m2 

Root  Chord 

4.5 

m 

Tip  Chord 

2.2 

m 

Aspect  Ratio 

17.5 

Wing  Incidence  Angle 

3.0° 

Horizontal  Tail  Span 

18.0 

m 

Horizontal  Root  Chord 

3.5 

m 

Horizontal  Tip  Chord 

2.45 

m 

Horizontal  Tail  Incidence  Angle 

-4.5° 

Vertical  Tail  Span 

4.0 

m 

Vertical  Root  Chord 

2.45 

m 

Vertical  Tip  Chord 

2.0 

m 

Wing/Horizontal  Tail  Airfoil 

NACA  4415 

Vertical  Tail  Airfoil 

NACA  0012 

Aileron  Location 

16.3  to  22.8 

m 

Aileron,  Elevator,  Rudder  Chord 

0.2  Ciocal 

Elevator  Span  Location 

1.8  to  9.0 

m 

Rudder  Span  Location 

0.8  to  3.2 

m 

Elements  per  wing 

9 

Elements  per  horizontal  tail 

5 

Elements  per  vertical  tail 

5 

Elements  in  fuselage 

10 

Total  Number  of  Elements 

48 

— 

Number  of  second-order  states 

192 

Number  of  first-order  states 

241 
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Table  10.  Control  and  aircraft  states  and  inertia  properties  of  the  very  flexible  aircraft  model 


Model  Parameters 

Property 

Value 

Units 

Light, 
no  thrust 

Heavy, 
no  thrust 

Light, 
thrust  for 
level  flight 

Heavy, 
thrust  for 
level  flight 

Elevator  deflec. 
angle,  Se 

4.51° 

-16.80° 

-6.89° 

-13.68° 

— 

Thrust  rqd., 

<5, 

0 

0 

3.21  •  104 

1.12  *  105 

N 

Aircraft  angle 
of  attack,  a 

0.80° 

7.62° 

1.93° 

7.30° 

Fuel  mass 

0 

32,000 

0 

32, 000 

kg 

Total  mass 

2.10-  104 

5.38  •  104 

2.10- 104 

5.38  •  104 

kg 

Fuel  fraction 

0.0 

59.5 

0.0 

59.5 

% 

TSS  * 

1  XX 

1.48-10° 

1.75-10° 

1.49-10° 

1.75-  10° 

kg  •  m2 

TSS 

yy 

8.20  •  105 

2.93  •  10° 

8.19  105 

2.93  •  10° 

kg  •  m2 

TSS 

1  ZZ 

2.27  •  10° 

4.46  •  10° 

2.26  •  10° 

4.47  •  10° 

kg  •  m2 

TSS 

xy 

0 

0 

0 

0 

kg  •  m2 

TSS 

JXZ 

0 

0 

0 

0 

kg  •  m2 

TSS 

*yz 

1.82  •  104 

9.20  •  104 

2.06  •  104 

9.00  •  104 

kg  •  m2 

xcm 

0 

0 

0 

0 

m 

Vcm 

3.13 

4.33  ■  10"3 

3.14 

5.64  •  10*3 

m 

zcm 

0.29 

0.79 

0.34 

0.77 

in 

*Note:  Isa  are  the  inertia  properties  in  a  deformed  steady  state  configuration 
Note:  All  aircraft  simulations  are  begun  at  sea  level  conditions  and 
65 m/s  level  flight 
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frame  and  in  the  y  direction,  such  that  a  positive  thrust  input,  6t  results  in  an  acceleration  in  the  positive  y 
direction. 

B.  Representative  Flexible  Aircraft  Trajectory  Tracking 

Flexible  aircraft  trajectory  tracking  studies  are  presented  here.  A  single  controller  is  presented  and  applied 
to  trajectory  control  for  climb  only,  bank  only,  and  simultaneous  climb  and  bank  all  at  a  given  (heavy) 
fuel  state.  The  controller  is  then  applied  at  an  alternate  (empty)  fuel  state  demonstrating  its  robustness  to 
significant  mass  and  inertia  changes.  Finally  a  representative  mission  profile  is  simulated. 

1.  Flexible  Aircraft  Controller  Design 

The  two  separate  inner-loop  controllers  for  a  flexible  aircraft  were  presented  in  Section  Ml.D.I.  For  the 
longitudinal  motion,  a  dynamic  inversion  approach  was  developed.  There  the  longitudinal  states  used  in 
development  of  the  controller  were  longitudinal  and  vertical  velocity  ( vbv  and  vbz)  and  pitch  rate,  t Jsr- 
The  corresponding  control  inputs  were  elevator,  8e ,  and  throttle,  6t.  Initially  all  three  longitudinal  states 
were  attempted  to  be  controlled  through  the  use  of  a  pseudo-inverse  of  the  resulting  control  effector  matrix 
function, ,  Eq.  7  *.  The  use  of  the  pseudo  inverse  resulted  in  a  stable  inner-loop  longitudinal  con¬ 
troller,  but  it  had  poor  tracking  performance.  A  significant  improvement  is  achieved  by  only  controlling  vbv 
and  ubx  •  The  resulting  set  of  error  states  are  modified  from  the  proposed  set,  Eq.  !)d  to 


Clong 


VBV  V  Bv  com 
UBX  -  UBrcom 


The  resulting  inner-loop  longitudinal  state  vector  is 


f  VBv 

Xiong  =  \  UJ bt  >  • 
\  &long  j 


(97) 


(98) 


A  standard  LQR  controller  is  applied  to  the  resulting  LTI  state  space  system  (see  Section  ■  for  details)  where 
the  weighting  matrices  are 


Qlong  =  diag[  1  1  102  104  ] 

(99) 

Rlong  =  diag  1  1  ]  • 

(100) 

For  the  lateral  inner-loop  controller,  the  error  states  are  also  modified  from  Eq.  1( 

H  to 

(101) 

^ By  OJ By  com  J 


This  was  done  for  two  reasons.  First,  due  to  the  addition  of  the  strain,  6,  and  strain  rate,  £,  states  to 
the  lateral  inner-loop  state  vector,  Eq.  h  i,  it  is  difficult  to  find  weighting  matrices  Q  and  R  which  provide 
solutions  to  the  ARE  of  the  LQR  technique.  Removal  of  the  yaw  rate  state  assisted  this  problem.  Second,  it 
was  found  that  commanding  the  states  of  lateral  velocity,  vbx,  and  roll  rate,  ujbv,  provided  stable  tracking 
of  the  removed  yaw  rate,  ljbz •  As  described  in  Section  1 11.0.2,  a  LQR  controller  is  applied  to  the  lateral 
motion  inner-loop  states.  The  weighting  matrices  are 


Qlat 

—  diag  ^  Qt  QvHt  Qu>Hv  Q^n,  Qelnt  ] 

(102) 

Rlat 

.2 

II 

(103) 

25  of  1.5 


American  Institute  of  Aeronautics  and  Astronautics 


where 


Q«  =  10-3, 

Qvh ,  =  0, 

=  10,  (104) 

Q^h,  =  o, 

Qeut  =  [  102  103  ]  • 

The  outer-loop  gains,  described  in  Section  111.1)..'),  are  loosely  designed  using  the  guidelines  of  Ziegler 


and  Nichols.44,45  The  final  chosen  values  are 

kq>  =  0.9;  =  0.225;  kd<t>  =  0.45,  (105) 

k<p  =  0.9;  kjQ  =  0.225;  kn4f  =  0.0225,  (106) 

A:7  =  0.5;  kh  =  0.5;  kJIl  =  0,  kdl  =  0.05;  (107) 

k^  =  0.5;  kjj  =  0.5;  fc//7  =  0.1.  (108) 


Finally  due  to  high-frequency  numerical  error  occurring  during  the  first  few  time  steps  of  the  simulations, 
a  third-order  low-pass  butterworth  filter56  was  applied  to  the  difference  equations  of  flight  path  angle,  7, 
and  Euler  bank  angle,  <F,  and  their  rates  given  in  Eq.  60.  The  filtered  equation  for  the  commanded  flight 
path  elevation  angle  difference  is 


X1  com  —  ^ fXlc.om  +  ^ flcom 

(7  com)/  —  CfXlcorn- 

For  the  third-order  2 -Hz  low-pass  butterworth  filter  used,  the  state  space  matrices  are 


Af 


Bf 

Cf 


-7T  0  0 

7T  -71  — 7T  , 

0  7T  0 

f 

7 r 

<  0  >, 

0 

[ooij- 


The  resulting  outer-loop  flight  path  and  Euler  bank  angle  commands  are  modified  from  Eq.  01  to 

A7com  =  MA7)/  +  kh  J"  (A7 )  dr  +  kill  as:  (A7)  dOj  dr  +  fcd7( A-y)/, 
Alcorn  =  fc7(A7)/  +  kiy(A-y)f  +  k,^  (A-y)  dr  +  knii  as:  (A7)  dr^  dr , 


A<I>com  =  A:<!>(A4>)/  +  fc/d>  /  (A4>)  dr  -f  kd,pA&f, 


f< 


A^con,  —  ^4>(A4>)/  +  kj<p(A$>)j  -+■ 


[  (A<J>)dr. 

Jo 


(109) 


(110) 


(111) 
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Figure  6.  Commanded  flight  path  angle  and  altitude,  heavy  weight  condition,  climb  command 


2.  Climb  Only 

Here  the  flexible  aircraft  is  commanded  to  make  a  100  m  climb  in  31  s  following  a  modified  cosine  input. 
Figure  <>  shows  that  both  the  altitude  and  flight  path  angle  have  excellent  tracking  to  the  commands.  The 
steady  altitude  error  is  less  than  0.5  m.  The  longitudinal  control  inputs  are  shown  in  Figure  7.  Other  than 
high-frequency  oscillations  due  to  numerical  initialization  error  seen  in  the  elevator  command,  Figure  7a, 
the  longitudinal  control  inputs  are  reasonable  and  reach  steady  state  values  upon  completion  of  the  climb. 
The  longitudinal  states  are  shown  in  Figure  >  and  0  along  with  the  desired  and  commanded  values.  Here 


time,  s  time,  s 

(a)  Elevator  (b)  Throttle 


Figure  7.  Longitudinal  flight  controls,  heavy  weight  condition,  climb  command 

excellent  tracking  is  seen  for  the  longitudinal  velocity,  vbv-  For  the  vertical  velocity  there  is  no  commanded 
value,  only  a  desired  trajectory  derived  from  the  desired  altitude  trajectory,  Figure  »»b,  and  the  nonlinear 
transformations  described  in  Section  lll.D.'t.  The  closed  loop  controller  provides  a  stable  vertical  velocity, 
as  seen  by  its  asymptotic  convergence  to  the  desired  steady  state  value.  Finally  the  aircraft  response,  inner 
loop  command,  and  desired  pitch  rate,  ljbt,  are  shown  in  Figure  T.The  effect  of  the  low-pass  butterworth 
filter  is  seen  by  examining  the  commanded  pitch  rate  and  the  simulated  pitch  rate  response.  Overall  the  use 
of  the  dynamic  inversion  on  the  longitudinal  inner- loop  dynamics  is  shown  to  work  very  well.  Additionally 
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Longitudinal  velocity, 


Figure  8.  B  reference  frame  velocities,  heavy  weight  condition,  climb  command 


Figure  9.  B  reference  pitch  rate  actual,  desired,  and  commanded  for  heavy  weight,  climb  command 
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this  example  shows  that  the  assumption  that  longitudinal  dynamics  may  be  decoupled  from  the  elastic  states 
was  valid  for  the  representative  aircraft. 

3.  Bank  Only 

Next  a  commanded  change  in  Euler  bank  angle  of  20°  while  maintaining  level  flight  is  simulated.  The  desired 
trajectories  and  aircraft  response  of  flight  path  angle,  altitude,  and  bank  angle  are  seen  in  Figure  10.  It  is 
immediately  apparent  that  the  longitudinal  and  lateral  coupling  is  excited  by  the  subsequent  change  in  flight 
path  angle  and  altitude.  However,  despite  an  almost  25  m  change  in  altitude,  the  altitude  is  controlled  via  the 
flight  path  angle.  Additionally  this  change  in  altitude  is  within  commercial  pilot  procedures,  and  subsequently 
would  not  create  a  problem  if  the  aircraft  were  flying  in  the  National  Airspace  system.  The  corresponding 


Figure  10.  Commanded  flight  path  angle,  altitude,  and  Euler  angle  for  heavy  weight  condition,  bank  command 


longitudinal  and  lateral  control  inputs  are  shown  in  Figures  I  I  and  12.  Steady  state  changes  are  seen  in  the 


Figure  11.  Longitudinal  flight  controls,  heavy  weight  condition,  bank  command 

elevator  and  throttle  from  an  initial  wings  level  values  of  —13.68°  and  1.12  •  10 RJV,  Table  0,  to  approximately 
-16.5°  and  1.25  •  10 5  TV.  Figure  12  shows  that  the  lateral  control  inputs  tend  to  steady  state  values  required 
to  maintain  a  level  turn  configuration.  The  controller  also  performed  well  in  smoothly  controlling  the  lateral 
motion  states  and  the  longitudinal  velocity  as  seen  in  Figures  13,  1 1,  and  15.  The  dynamic  inversion  inner- 
loop  controller  again  shows  excellent  tracking  between  commanded  and  simulated  response,  Figure  1 3b  and 
15a.  The  oscillations  of  the  longitudinal  motion  about  desired  trajectories  (Figures  13a,  b,  and  15a)  are 
attributed  to  the  difficulty  in  tuning  the  nonlinear  outer-loop  controller  when  a  lateral  motion  (bank  angle 
here)  is  commanded.  However  despite  this  difficulty,  the  controller  does  exhibit  stable  tracking. 
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Longitudinal  velocity, 


Figure  12.  Lateral  flight  controls,  heavy  weight  condition,  bank  command 


(a)  Longitudinal  velocity 


Figure  13.  B  reference  frame  velocities,  heavy  weight  condition,  bank  command 
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Figure  14.  B  reference  frame  lateral  velocity,  heavy  weight  condition,  bank  command 


time,  s  time,  s 


(a)  Pitch  Rate  (b)  Roll  Rate 


(c)  Yaw  Rate 


Figure  15.  B  reference  frame  angular  velocities,  heavy  weight  condition,  bank  command 


31  of  13 


American  Institute  of  Aeronautics  and  Astronautics 


4-  Bank  arid  Climb 


Finally  for  the  heavy  weight  fuel  condition,  a  combined  climb  and  turn  trajectory  is  commanded.  The 
flight  path  angle,  altitude,  and  bank  angle  change  trajectories  are  seen  in  Figure  16.  The  corresponding 


Figure  16.  Commanded  flight  path  angle,  altitude,  and  Euler  angle  for  heavy  weight  condition,  bank/climb 
command 


longitudinal  and  lateral  control  inputs  are  shown  in  Figures  1 7  and  I  n.  The  B  reference  frame  linear  and 


Figure  17.  Longitudinal  flight  controls,  heavy  weight  condition,  bank/climb  command 

angular  velocities  are  shown  in  Figures  Ml,  20,  and  21.  Similar  difficulties  to  the  bank  only  command 
are  seen  in  this  simulation.  However,  due  to  the  overall  good  control  on  total  velocity,  Eq.  27,  Figure  10 
(between  64.8  to  65.2™),  and  acceptable  control  on  the  flight  path  angle,  Figure  10,  there  is  less  than  a  5  rn 
steady  state  error  in  altitude.  Additionally,  the  desired  Euler  bank  angle  is  achieved  at  steady  state. 

5.  Bank  and  Climb  at  Empty  Fuel  State 

To  investigate  the  robustness  of  the  proposed  controller  architecture,  the  same  bank  and  climb  trajectory 
of  the  previous  section  is  repeated  here,  but  at  an  empty  fuel  condition.  Figure  22  shows  the  desired  and 
actual  response  in  flight  path  angle,  altitude,  and  bank  angle.  While  the  bank  angle  response  is  similar  to 
heavy  fuel  condition,  Figure  16c,  there  are  greater  excursions  in  the  flight  path  angle  and  altitude.  These 
excursions  are  attributed  to  the  outer-loop  controller.  The  longitudinal  and  lateral  control  inputs  are  shown 
in  Figures  22  and  2  1.  The  B  reference  frame  linear  and  angular  velocity  trajectories  (actual,  commanded, 
desired)  are  shown  in  Figures  25,  26,  and  27.  An  examination  of  the  vertical  velocity,  Figure  27b,  highlights 
a  problem  with  applying  the  controller  designed  for  the  heavy  weight  to  the  light  weight  condition.  In  the 
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Longitudinal  velocity,  v8/  m/s  Aileron  Deflection  deg 


Figure  18.  Lateral  flight  controls,  heavy  weight  condition,  bank/climb  command 


Figure  19.  B  reference  frame  velocities,  heavy  weight  condition,  bank/climb  command 
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4 


Figure  20.  B  reference  frame  lateral  velocity,  heavy  weight  condition,  bank/climb  command 


Figure  21.  B  reference  frame  angular  velocities,  heavy  weight  condition,  bank/climb  command 


Figure  22.  Commanded  flight  path  angle,  altitude,  and  Euler  angle  for  light  weight  condition,  bank/climb 
command 
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Deflection,  deg  Elevator  Deflection,  deg 


Figure  23.  Longitudinal  flight  controls,  light  weight  condition,  bank/climb  command 


Figure  24.  Lateral  flight  controls,  light  weight  condition,  bank/climb  command 


35  of  15 


American  Institute  of  Aeronautics  and  Astronautics 


outer-loop  nonlinear  transformation,  Section  III.D.3,  the  body  angle-of-attack,  a,  is  a  derived  and  required 
parameter,  based  upon  a  steady  state  wings-level  zero  angle-of-attack,  q0,  given  in  Eq.  Is.  The  heavy  weight 
controller  was  designed  using  a  trim  body  body  angle  of  attack,  q()  of  7.30°  and  presented  in  Table  Id,  and 
applied  to  the  light  weight  condition  here.  Note,  for  a  light  weight  condition,  the  trim  body  angle  of  attack, 
Qo>  is  1.93°.  Despite  the  resulting  large  discrepancy  in  q0  and  vertical  velocities  between  the  light  and  heavy 
weight  conditions,  the  controller  performs  adequately  for  the  significant  change  in  mass,  Table  Id.  This 
is  because  vertical  velocity  is  not  actively  controlled  and  the  other  B  reference  frame  velocities  are  not  as 
significantly  affected  by  a(). 


£ 

i 

a> 

> 


time,  s  time,  s 


(a)  Longitudinal  velocity 


(b)  Vertical  velocity 


Figure  25.  B  reference  frame  velocities,  light  weight  condition,  bank/climb  command 


Figure  26.  B  reference  frame  lateral  velocity,  heavy  weight  condition,  bank/climb  command 


6.  Mission  Profile  Segment  Case 

A  final  representative  mission  profile  segment  simulation  is  presented  here.  The  aircraft  is  commanded 
to  perform  the  following  sequence  of  maneuvers:  climb  and  turn,  roll  out  of  turn  and  continue  climbing, 
level  off,  level  turn,  descent  and  turn,  roll  out  of  turn  and  continue  descending,  and  level  off  as  seen  in 
Figure  -  V  The  altitude  and  associated  flight  path  angle  along  with  the  Euler  bank  angle  desired  and  actual 
trajectories  are  seen  in  Figures  29  and  30.  Altitude  tracking  presents  a  maximum  overshoot  of  46  meters 
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Figure  27.  B  reference  frame  angular  velocities,  light  weight  condition,  bank/climb  command 


following  the  commanded  climb.  This  overshoot  occurs  while  the  bank  angle  is  settling  to  its  commanded 
zero  degree  angle.  The  overshoot  in  altitude  is  acceptable  given  that  this  class  of  aircraft  is  expected  to 
fly  in  the  National  Airspace  System.  Typical  piloted  aircraft  flying  under  instrument  flight  rules  (IFR)  are 
required  to  maintain  less  than  about  30  meter  steady  state  altitude  deviations  and  are  allowed  up  to  60  meter 
temporary  deviations.  The  corresponding  longitudinal  and  lateral  control  inputs  are  seen  in  Figures  3 1  and 
32.  During  the  second  portion  of  the  commanded  climb,  an  increase  in  throttle  and  decrease  in  elevator 
are  seen  (between  50  and  75  seconds)  corresponding  to  the  greater  commanded  climb  rate  (6.67 m/s  at  25 
seconds  and  13.33m/s  at  65  seconds).  There  is  a  steady  state  altitude  error  of  less  than  2  meters  at  250 
seconds  and  4  meters  at  400  seconds.  Bank  angle  overshoots  are  less  than  7.2°  with  excellent  steady  state 
tracking.  The  B  reference  frame  desired,  commanded,  and  actual  velocities  are  seen  in  Figures  33,  3  1, 


Figure  28.  Representative  mission  profile  segment 

35,  and  36.  As  with  previous  results  there  is  excellent  tracking  of  the  commanded  longitudinal  velocity, 
Figure  33b,  and  pitch  rate,  Figure  35a.  Again  tracking  error  between  the  desired  and  commanded  values 
in  the  longitudinal  motion  are  attributed  to  outer-loop  controller  architecture.  For  the  lateral  motion  the 
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levator  Deflection,  dpg  ,  Flight  Path  Angle,  y,  deg 


Figure  29.  Commanded  flight  path  angle  and  altitude  for  heavy  weight  condition 


Figure  30.  Commanded  Euler  bank  angle  for  heavy  weight  condition 


Figure  31.  Longitudinal  flight  controls,  heavy  weight  condition 
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Figure  32.  Lateral  flight  controls,  heavy  weight  condition 


sideslip  velocity,  Figure  .Ida,  is  seen  to  be  below  0.04m/s  or  0.05%  of  the  total  velocity,  V.  Additionally  the 
roll  rate  (Figure  d.">b)  has  acceptable  tracking  and  the  yaw  rate  (Figure  .T»c),  while  not  directly  controlled, 
has  relatively  good  and  stable  tracking  of  the  desired  trajectory. 


(b)  Vertical  velocity 


Figure  33.  B  reference  frame  velocities,  heavy  weight  condition 


V.  Concluding  Remarks 

The  control  architecture  was  shown  numerically  to  track  altitude  and  bank  angle  changes  in  the  presence 
of  smooth  air,  full  elastic  state  feedback,  and  perfect  sensors.  The  purpose  of  these  studies  was  not  to  design 
a  perfect  controller,  but  rather  demonstrate  the  viability  of  the  proposed  method.  In  such,  there  is  more 
tuning  which  could  be  accomplished  for  tighter  trajectory  control.  The  inner-  and  outer-loop  architecture 
demonstrated  an  effective  means  for  accomplishing  trajectory  control  of  very  flexible  aircraft.  While  more 
complex,  the  proposed  architecture  shows  to  be  much  easier  to  achieve  stable  tracking  as  compared  to  more 
traditional  methods.  Furthermore,  the  separation  of  longitudinal  and  lateral  motion  control  handled  by 
separate  controllers  proved  to  be  very  effective.  Finally,  while  not  presented  here,  the  general  trends  of 
increased  integral  gain  destabilizing  system  response  and  increased  derivative  gain  stabilizing  and  slowing 
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Pitch  Rate,  to  e ,  1 0-2  rad/s  Lateral  Velocity,  ve  ,  10-2  m/s 


time,  s  time,  s 

(a)  Lateral  velocity  (b)  Lateral  velocity  as  a  percentage  of  total  velocity 

Figure  34.  B  reference  frame  lateral  velocity,  heavy  weight  condition 


Figure  35.  B  reference  frame  angular  velocities,  heavy  weight  condition 


Figure  36.  Yaw  rate,  B  reference  frame  angular  velocity,  heavy  weight  condition 


40  of  M 


American  Institute  of  Aeronautics  and  Astronautics 


system  response  were  seen  in  the  development  of  the  outer-loop  nonlinear  gains. 
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RQ-4A/B  GLOBAL  HAWK  HIGH  ALTITUDE,  LONG 
ENDURANCE,  UNMANNED  RECONNAISSANCE  AIRCRAFT, 
USA 


RQ-4A  Global  Hawk  is  a  high  altitude,  long  endurance  unmanned  aerial 
reconnaissance  system  which  provides  military  field  commanders  with  high 
resolution,  near-real-time  imagery  of  large  geographic  areas. 
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The  program  is  funded  by  the  Defense  Airborne 
Reconnaissance  Office  (DARO)  and  managed  by 
the  Defense  Advanced  Research  Projects  Agency 
(DARPA)  and  the  US  Air  Force. 


reconnaissance 

system." 


Northrop  Grumman  Corporation,  Ryan 
Aeronautical  Centre  is  the  prime  contractor  and  the 
principal  suppliers  include  Raytheon  Systems 
(sensors),  Rolls-Royce  North  America  (turbofan 
engine),  Boeing  North  American  (carbon  fibre  wing)  and  L3  Communications 
(communications  system).  The  Global  Hawk  air  vehicles  are  built  at  the 
Northrop  Grumman  (formerly  Teledyne  Ryan)  Aeronautical  facility  in  San 
Diego. 
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GLOBAL  HAWK  DEVELOPMENT 

In  March  2001,  the  US  Department  of  Defense  awarded  Northrop  Grumman  a 
contract  for  the  Engineering  and  Manufacturing  Development  (EMD)  phase  of 
the  program  which  concluded  in  February  2003  with  the  final  delivery  of  the 
seventh  pre-production  (Block  0)  vehicle. 
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In  June  2001  a  contract  was  placed  to  begin  Low  Rate  Initial  Production  (LRIP) 
for  two  production  air  vehicles  and  the  mission  control  element  of  the  system's 
ground  station,  to  be  completed  by  December  2003. 
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The  first  production  vehicle  (Block  10)  rolled  out  in  August  2003.  A  further  LRIP 
contract  for  four  vehicles  was  placed  in  February  2003  and  a  third  in  October 
2004  for  two  vehicles.  Block  10  deliveries  were  completed  in  June  2006. 


The  US  Navy  had  two  air  vehicles  delivered  in  2005. 
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Northrop  Grumman  is  developing  the  next- 
aeneration.  RQ-4B.  which  has  a  50%  Davload 
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